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1.  Summary 


It  is  well-known  that  low  frequency  electromagnetic  (EM)  signals  are  heavily  attenuated  in 
a  medium  with  dense  electron  population.  All  signals  below  the  plasma  frequency  of  the 
medium  get  cut  off.  If  we  create  in  this  medium  a  small  population  of  relatively  hot 
electrons  the  composite  medium  then  supports  low-frequency  electrostatic  oscillations 
known  as  electron  acoustic  waves  (EAW)  [e.g.,  Gary  and  Tokar,  Phys.  Fluids,  vol.  28,  p 
2439].  The  dispersion  relation  of  this  composite  medium  shows  that  it  supports  EAW  in 
the  frequency  band  where  EM  signals  are  cut  off.  Thus  it  is  possible,  in  principle,  to 
employ  EAW  to  transmit  signals  across  an  overdense  plasma  medium.  Our  primary 
interest  in  this  report  is  to  study  the  radiation  characteristics  of  a  source  current  distribution 
embedded  in  a  half-space  of  our  composite  medium.  To  enable  this,  we  derive  the  Green's 
functions  for  our  problem  and  hence  study  the  radiation  characteristics  of  antennas.  When 
the  source  signal  frequency  is  below  the  plasma  frequency,  only  EAW  exist  in  the 
composite  medium,  while  only  EM  waves  can  exist  in  the  free  space  above.  We  find  that 
the  far-zone  radiation  fields  of  any  current  distribution  consist  only  of  9-polarized  waves. 
Explicit  expressions  for  the  radiated  fields  are  obtained  for  horizontally-  and  vertically- 
polarized  Hertzian  dipoles  embedded  in  our  composite  medium.  We  hence  find  that  in  both 
cases  the  radiation  patterns  are  skewed  towards  the  horizon.  In  particular,  we  find  that  the 
radiation  pattern  of  a  horizontal  dipole  has  two  lobes  as  opposed  to  one  in  the  underdense 
case. 
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2.  Introduction 


When  a  space  vehicle  reenters  the  atmosphere  it  encounters  a  bow  shock.  The  kinetic  energy  of 
the  vehicle  is  translated  into  high  heat  and  pressure  in  the  shock  layer.  This  initiates  a  variety  of 
thermo-chemical  processes  [1],  [2]  and  a  dense  plasma  sheath  surrounds  the  vehicle.  It  is  well- 
known  that  sensor  performance  aboard  hypersonic  vehicles  is  degraded  to  various  degrees  because 
of  the  dense  plasma.  In  our  study  we  have  found  that  there  will  be  periods  in  the  course  of  the 
trajectory  where  the  electron  density  of  the  plasma  sheath  is  so  high  that  the  microwave  signals  are 
completely  cut  off.  This  is  the  famous  "communication  blackout"  encountered  during  the  Apollo 
missions  in  the  1960s.  Indeed  issues  with  communication  with  reentry  vehicles  were  studied  even 
before  Apollo  missions  [3].  Such  communication  blackout  problems  are  also  encountered  in  many 
other  space  exploration  missions  [4],  Some  of  early  schemes  that  were  proposed  are:  use  of  electron 
beams  [5],  high  frequency  laser  [6],  strong  static  magnetic  fields  [5]-[8],  and  alteration  of  chemical 
composition  [9]-[13]. 

Two  articles  that  summarize  early  ideas  are  those  of  Rawhouser  [5]  and  Rybak  [14].  More 
recent  ideas  put  forward  are:  three-wave  scattering  process  [15],  [16],  Hall  effect  drift  [17], 
electrostatic  manipulation  [  1 7]-[  1 9],  and  resonant  transmission  [20],  [21].  These  schemes  exploit 
the  wave  phenomenology  in  the  plasma  sheath  to  transmit  signals  across  the  sheath.  There  have  also 
been  suggestions  on  mechanical  structures  that  can  avoid  the  blackout  condition.  One  idea  is  that  of 
a  remote  antenna  assembly  [22]  just  outside  the  potential  dense  plasma  region.  Also,  the  choice  of 
shape  of  the  hypersonic  vehicle  has  an  important  influence  on  the  character  of  the  flow  field  [23]. 
Hence  it  is  possible  to  choose  the  shape  of  the  vehicle  such  that  blackout  situations  are  minimized. 
All  these  schemes  have  their  share  of  advantages  and  disadvantages.  None  have  been  implemented 
successfully  in  any  hypersonic  vehicles  thus  far.  Indeed  there  is  a  long  way  to  go  before  any  of 
these  schemes  can  be  implemented  in  operational  hypersonic  vehicles  because  there  remain  many 
theoretical  and  practical  issues  to  be  overcome.  Some  of  the  difficulties  one  faces  are:  (a)  size, 
weight  constraint,  (b)  engineering  issues  of  implementation,  (c)  power  resources  constraint. 
Furthermore,  the  proposed  concepts  have  been  illustrated  only  for  idealized  models;  the  feasibility 
and  success  of  such  concepts  are  questionable  under  non-ideal  conditions  encountered  in  practice. 
Hence  the  issue  of  communication  blackout  is  as  "topical"  as  it  was  during  Apollo  missions  [24], 
[25].  Although  there  are  now  other  means  to  communicate  (via  satellite)  the  sensor  needs  for 
modern  hypersonic  vehicles  are  now  much  greater  and  more  critical.  One  concept  that  we  are 
currently  investigating  is  to  employ  electron  acoustic  waves  (EAW)  for  carrying  low-frequency 
microwave  signals  across  the  overdense  plasma  sheath  [26].  The  idea  is  to  create  a  small  population 
of  high  temperature  electrons  into  the  plasma  sheath  so  that  we  have  a  two-temperature  electron 
population.  One  way  to  create  this  composite  medium  is  by  injecting  hot  electrons  by  an  electron 
gun  as  it  was  done  during  CHARGE-2B  experiments  [27]-[29] .  Although  EAW  have  not  been 
popularly  studied  in  the  literature  there  are  a  few  interesting  papers  [28],  [30]  where  they  play 
important  roles  in  the  context  of  ionosphere  and  space  experiments.  Indeed  there  are  several 
instances  where  a  two-temperature  electron  population  occurs  in  nature.  Some  examples  are:  earth's 
bow  shock  [31],  [32],  undisturbed  solar  wind  [35],  interplanetary  shocks  [36],  In  these  contexts  the 
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concept  of  EAW  are  sometimes  employed  to  offer  theoretical  explanation  to  certain 
phenomenology. 

Note  that  it  is  not  mandatory  to  have  a  two-temperature  electron  population  to  excite  EAW.  One 
may  as  well  have  a  single  electron  population  along  with  a  resonant  electron  beam  to  excite  EAW 
[37].  Although  we  have  not  explored  this  path  in  detail  we  feel  that  the  domain  of  applicability  of 
our  scheme  will  be  much  larger  than  that  using  single  electron  population.  These  remarks  are  based 
on  the  assumption  that  the  electrons  have  Maxwellian  distribution.  Furthermore,  there  are  other 
ways  (see  for  example,  [37],  [38]),  besides  that  using  a  resonant  electron  beam,  for  creating  a 
plateau  near  the  region  of  phase  velocity  of  EAW.  However,  not  all  such  schemes  are  suitable  for 
use  in  hypersonic  vehicles.  Since  we  have  more  familiarity  and  confidence  in  the  beam  concept  for 
our  application,  our  study  is  restricted  to  this  idea. 

In  this  report  we  show  that  the  two-temperature  electron  population  supports  electron  acoustic 
waves  [3  9] -[41]  in  the  frequency  domain  where  the  low-frequency  electromagnetic  waves  are  cut 
off.  In  fact  the  dispersion  relation  shows  that  the  propagation  characteristics  of  EM  waves  and  EAW 
complement  each  other.  Thus  it  is  a  promising  idea  to  employ  EAW  to  carry  low-frequency 
microwave  signals  across  an  over  dense  plasma  slab.  However,  there  are  issues:  one  of  them  is  the 
phenomenon  of  resonant  damping.  A  kinetic  approach  is  taken  to  understand  the  damping 
characteristics  of  EAW.  We  also  investigate  the  radiation  characteristics  of  Hertzian  dipole  antennas 
embedded  in  the  two-temperature  composite  medium.  In  this  connection  we  derive  the  Green's 
function  for  our  overdense  plasma  layer.  We  hence  find  that  the  radiation  characteristics  of  antennas 
in  an  overdense  two-temperature  electron  population  are  quite  different  from  those  in  free  space:  (a) 
only  6  -polarized  waves  are  radiated,  (b)  radiation  patter  is  skewed  towards  the  horizon  and  (c)  the 
radiation  pattern  of  a  horizontal  dipole  has  two  lobes  as  opposed  to  one  in  the  underdense  case. 
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3.0  Methods,  Assumptions,  and  Procedures 

3.1  Electromagnetic  Waves 


First  we  study  the  characteristics  of  electromagnetic  waves  (EMW)  in  our  problem.  Maxwell's 
equations  are  given  as 

VxH  +  ^E-— j  =  — J  (1) 

c  c  J  c  v  7 

VxE  =  fH  (2) 

Here  J  is  the  source  current  exciting  the  system;  j  is  the  current  distribution  due  to  the  plasma 
sheath  given  as 

j  =  -eno\  (3) 

where  e  is  the  electron  charge  and  no  is  the  unperturbed  electron  charge  density.  The  velocity  v  of 
the  electrons  is  governed  by  the  following  linearized  equation  of  motion: 

a,v  +  ^  =  -^E  (4) 

t  mn0  m  x  7 

where  T  is  the  electron  temperature  and  m  is  the  electron  mass.  Note  that  we  have  neglected  here 
the  electron-neutral  collisions 1 .  The  density  and  velocity  of  the  electrons  are  connected  by  the 
continuity  equation 

dtn  +  V  •  ( na\ )  =  0 

On  taking  the  curl  of  (2)  and  inserting  (1)  in  it  we  obtain 


VxVxE-£2E-^^v  =  ^J 

c  c 

where  k  —  co/c  is  the  free  space  wave  number.  From  the  equation  of  motion,  we  have 

v  =  -r^E 

icom 

Substituting  (7)  in  (6)  we  obtain  the  following  equation  for  the  electromagnetic  waves 

V  x  V  x  E  -  k2E  =  ^  J 

e  C 

where 

k2e=k2[l- 

and  <x>p  is  the  plasma  frequency  given  as 

col 


4  7ie2nn 


(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


1 


Collisions  will  not  affect  the  fundamental  character  of  EMW  or  the  electron  acoustic  waves.  They  will  only  add  attenuation  to  them. 
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When  co>  cop,  then  k2  is  less  than  k2  and  is  positive. This  means  that  the  medium  will  support 
propagating  electromagnetic  waves.  On  the  other  hand  when  co  <  cop  ,  k2e  is  negative  and  the 
electromagnetic  waves  here  become  evanescent  .  The  dispersion  relation  of  the  EMW  is  given  as 

co2  =  orp  +  k2c2  (11) 

Thus  we  see  that  cop  is  the  cutoff  frequency.  All  signals  with  frequency  components  below  cop  are 
cut  off.  Notice  that  cop  is  proportional  ,  and  hence,  the  cutoff  frequency  increases  with 
increasing  density.  The  dispersion  relation  is  graphically  displayed  in  Figure  1. 


Fig.  1.  Dispersion  relation  of  electromagnetic  waves  for  various  electron  densities 

We  have  plotted  frequency  versus  wave  number  for  various  electron  densities  (cirT3 * ).  We  have 
considered  electron  densities  in  the  range  106  to  1016 .  The  range  of  parameters  that  describe  the 
hypersonic  environment  is  very  large.  Several  factors  such  as  velocity,  altitude,  trajectory,  shape 
and  size  of  the  vehicle,  and  heat  shield  materials  have  strong  influence  on  the  environment.  Notice 
that  for  n  —  1 010 ,  GPS  signals  are  cut  off;  at  n  —  1 012 ,  communications,  telemetry,  and  radar  signals 
are  cut  off  .  Therefore  the  problem  that  we  may  encounter  is  quite  serious.  Although  several 
remedies  have  been  suggested,  all  have  their  share  of  limitations  and  shortcomings.  For  instance, 
there  are  engineering  issues,  size-weight  constraints,  and  system  integration  issues.  Also,  most 
mitigation  ideas  have  only  been  proposed  and  illustrated  by  assuming  simplified  ideal  conditions. 
We  are  currently  exploring  the  idea  of  employing  the  electron  acoustic  waves  in  a  two-temperature 
electron  population.  One  can  inject  a  small  population  of  hot  electrons  such  that  we  have  a 
composite  medium  consisting  of  cold  and  hot  electrons.  We  find  that  such  a  composite  medium  can 


2  This  situation  is  called  the  overdense  condition.  The  other  situation  when  co>  co  is  called  the  underdense  condition. 

P 

^  The  interest  in  this  report  is  the  domain  when  the  thickness  of  the  plasma  sheath  is  much  larger  than  the  skin  depth  of  the  EMW. 
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support  stable  electron  acoustic  waves.  We  see  that  the  propagation  characteristics  of  EAW  in  our 
medium  complement  those  of  EMW  in  the  sense  that  EAW  propagate  in  the  frequency  band  where 
EMW  are  cut  off.  Hence  it  is  a  good  idea  to  employ  EAW  to  carry  low-frequency  EMW  across  an 
overdense  plasma  sheath. 


3.2  Electron  Acoustic  Waves 


We  employ  a  hydrodynamic  description  of  the  two-temperature  electron  population  problem. 

The  governing  equations  are  the  same  as  in  (1),  (2),  (3),  and  (5).  But  now  the  current  due  to  the 
charge  distribution  has  two  components:  one  due  to  hot  electrons  and  the  other  due  to  cold 
electrons.  Hence 

j  =  —enocvc  —  enohvh  (12) 

where  subscripts  h  and  v  stand  for  hot  and  cold  electrons.  There  are  thus  two  linearized  equations 
of  motion  (one  for  hot  and  one  for  cold  electrons): 

—  =  -- £-E  a  =  \h,c\  (13) 

t  a  mnoc ■  m  l  ’  )  v  7 

Also,  there  are  two  continuity  equations  given  as 

dtna+V-(noa  vj  =  0  (14) 

We  apply  the  divergence  and  time-derivative  operator  on  the  first  Maxwell's  equation  to  get 

co2V  ■  E  +  4 nedt  {noV  ■  \c  +  nohV  ■  \h }  =  -4 moN  J  (15) 

We  have  assumed  that  the  unperturbed  charge  densities  noc  and  noh  are  constants  in  our  problem. 
For  the  case  of  cold  electrons,  we  assume  that  Tc  is  small  and  neglect  the  tenn  involving  T  in  the 
equation  of  motion.  Thus 

<9  v  □  E  (16) 

t  c  m  v  7 

Further  note  that  EAW  are  longitudinal  and  hence  irrotational.  Hence  we  express  E  as  the  gradient 
of  the  scalar  potential  y/  .  Thus 

4 7venocdtV  ■  \c  =  -co2cV2y/  (17) 

where  coc  is  the  plasma  frequency  of  the  cold  electrons.  For  the  hot  electrons  we  use  the  continuity 
equation  to  express  V  •  \h  and  use  the  following  relation  obtained  from  the  linearized  equation  of 
motion: 


Thus 


ii, 

mn„ 


-Vnh  □  -- 


E 


4nenoh8y  -\h  =-^y/ 

Hi 


where  rh  is  the  Debye  radius  of  the  hot  electrons  defined  as 


(18) 

(19) 

(20) 
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where  Vh  and  coh  are  the  thermal  speed  and  plasma  frequency  of  the  hot  electrons.  On  substituting 
(17)  and  (19)  in  (15)  we  obtain  the  following  equation  for  EAW: 

VV  +  £>  =  4^V-  J  (21) 

'  a  1  0)c  —co 

with 


k 2 


a 


(22) 


where  ka  is  the  wave  number  of  EAW  in  our  composite  medium.  Thus  the  dispersion  relation  for 
EAW  is  given  as 


co2  = 


\+k 


-v 


(23) 


Notice  that  in  contrast  to  EMW  the  plasma  frequency  is  now  the  upper  cutoff  frequency.  In  other 
words,  the  medium  supports  EAW  only  with  frequencies  less  than  the  plasma  frequency  of  the  cold 
electron  population.  The  dispersion  relation  of  EAW  is  graphically  displayed  in  Figure  2.  We  have 
chosen  the  hot  electron  density  to  be  5%  of  the  total  electron  population  and  its  temperature  to  be 
20  eV.  Since  the  cold  electron  temperature  is  taken  as  0.2  eV  we  have  a  temperature  contrast  ratio 
of  100. 


Figure  2.  Dispersion  relation  of  electron  acoustic  waves  for  various  total  electron  densities 

It  is  instructive  to  compare  Figures  1  and  2.  Both  graphs  have  the  same  set  of  electron 
densities.  Notice  that  for  EMW  the  plasma  frequency  is  the  lower  bound  of  its  pass  band.  In 
contrast,  for  EAW  the  plasma  frequency  is  the  upper  bound  of  its  pass  band.  In  that  sense  they  are 
exact  complements.  However,  there  are  differences  in  the  manner  in  which  they  approach  cutoff.  In 

the  case  of  EMW  co/k  — »  c  for  coU  cope .  In  the  case  of  EAW,  co/k  — »  Vh  - \jnjnh 

for  coU  cope. 
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Two  basic  quantities  characterizing  waves  are  phase  velocity  ( v  )  and  group  velocity  ( vg  ), 
which  are  defined  as 


For  the  case  of  EMW, 


=f 


v  =  — 

P  v 


v  = 

Vg  9k 


v=  1- 


CO  >  CO, 


Since  v  <  1 ,  the  phase  velocity  of  the  EMW  is  larger  than  c  However  the  group  velocity 
is  vc  and  hence  is  less  than  c  . 


(24) 

(25) 


For  the  case  of  EAW, 


vg=vpvS 


5  = 


1 +k 


2  2 
h 


(26) 


EMW  may  be  considered  as  constant  velocity  waves  (as  is  evident  for  large  k  )  while  EAW  are 
constant  frequency  waves.  For  large  wave  numbers  the  phase  velocity  of  EMW  approaches  the 
constant  value  c  .  In  contrast  the  phase  velocity  of  EAW  decreases  with  increase  in  k  .  Thus  the 
group  velocity  of  EAW  vanishes  for  large  k  .  However,  we  will  see  that  EAW  are  weakly  damped 
only  for  krh  □  1  and,  hence,  we  will  operate  in  the  domain  where  the  group  velocity  is  nonzero  so 
that  we  can  carry  low-frequency  signals  across  the  overdense  plasma  sheath  using  EAW. 
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3.3  Green’s  Function  of  Electron  Acoustic  Waves 


In  order  to  understand  the  radiation  characteristics  of  EAW  we  now  derive  their  Green's 
function.  The  governing  equation  is  (21).  The  conservation  equation  for  source  current  there  is 
given  as 

V  •  J  +  dtq  =  0  (27) 

where  q  is  the  charge  density  of  the  source.  Inserting  this  relation  in  (21),  we  obtain 

V2¥  +  k2a¥  =  -47zqk2ar2h  (28) 
Employing  the  Fourier  transform  we  obtain  the  following  solution: 

¥(k,6))  =  -^47rq  (29) 

Thus 

¥(r,a)  =  j^^1^47rqe‘krdk  (30) 

To  evaluate  the  integral  we  introduce  polar  coordinates  k,0,(f)  as  k  .  =  /f  sin  Ocos(f) . 

kv  -  k  sin  $sin  <j)  ,  k_  -  k  cos  0  and  assume,  without  loss  of  generality,  that  r  =  R  z .  Thus  we  have 


¥{r, co)  =  ^ J ^  j e,kRcos0  sin OdOdk 
0  *  0 


(31) 


Since  EAW  undergo  resonant  damping  (as  will  be  seen  in  the  next  section)  ka  is  complex.  The 
integrand  thus  has  two  poles,  at  k  -  ±ka  as  shown  in  Figure  3.  Noting  that  e'kR  =  e'kRe~k  R  ( k'  and 


k"  are  the  real  and  imaginary  parts  of  k  ),  we  choose  the  contour  of  integration  along  the  semi¬ 
circle  in  the  upper  half  plane  as  shown  in  the  figure.  Since  the  integrand  vanishes  along  the  semi¬ 
circular  path  the  integral  in  (3 1)  becomes 


I 


k~elkR  7/  .  ik  R 

dk  =  me  ° 


k-r 


(32) 


where  C  denotes  the  contour  shown  in  Figure  3.  Substituting  this  in  (3 1),  we  obtain 


¥{x,co) 


_  ikaR 

R  ^ 


(33) 


From  this  we  obtain  the  space-time  potential  of  EAW  as 


¥(r,t)  =  i V  j  v{r,CD)emd(o 


(34) 
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The  propagation  constant  of  EAW  is  given  as 

k  c  =  r, 

a  c,,  5  a  h 


^Rr- 


0)~ 


(35) 


where  ca  is  the  phase  velocity  of  EAW  which,  as  shown  above,  varies  with  frequency. 

This  dispersive  character  makes  the  evaluation  of  (34)  rather  difficult.  However,  when  co  □  coc ,  the 
phase  velocity  is  approximately  constant  given  as 


r  —V  — 

Ca  V  h'>  - 


(36) 


(37) 


Now  we  may  proceed  to  evaluate  (34)  as 

oo 

^(r>0=  2^1  J  j  CO2  QKp{ia>(f  -t^d(o 

—00 

—*$*•(*- ') 

where  8 "  is  the  second  derivative  of  the  delta  function.Thus  the  space-time  Green's  function  for 
EAW  is  written  as 

=  (38) 

where  U  is  the  Heaviside  unit  step-function.  Note  that  this  form  of  Green's  function  for  the  EAW  is 
valid  only  under  the  approximation  that  the  phase  velocity  is  a  constant.  As  seen  in  Figure  2  the 
phase  velocity  is  indeed  a  constant  when  coU  coc .  Therefore,  if  the  signal  bandwidth  falls  well 

below  0J„  the  above  Green's  function  is  accurate. 
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3.4  Kinetic  Approach  to  Electron  Acoustic  Waves 


Our  study  of  EAW  in  the  last  section  was  based  on  a  hydrodynamic  description  of  plasma. 
Although  the  dispersion  relation  is  correct  it  does  not  account  for  the  phenomenon  of  resonant 
damping4  which  has  a  major  impact  on  the  ability  of  EAW  to  carry  low-frequency  microwave 
signals.  To  assess  this  impact  we  adopt  a  kinetic  approach.  The  fundamental  quantity  here  is  the 
velocity  distribution  function  which  is  governed  by  the  following  equation. 

dtf  +  v  •  V/  =  ^-E  •  Vv/  (39) 


We  have  one  such  Vlasov  equation  for  the  hot  and  cold  electron  populations.  Because  of  its  slow 
velocity,  ion  population  does  not  play  any  role  in  the  estimate  of  resonant  damping.  We  ignore 
collisions  and  use  a  linearized  Vlasov  equation  along  with  the  Poisson  equation.  Thus  we  arrive  at 
the  following  dispersion  equation 


1+7M1+ CZ(C)}  +  rfrfl  +  CZ(C)}  =  0 


(40) 


where  Z  is  the  plasma  dispersion  function  defined  as 

z(Cc)  = 

dy 

(41) 

and 

ii 

^ s 

a  =  (c,  /?} 

(42) 

EAW  exist  in  the  regime  where  the  phase  velocity  lies  well  within  the  interval  between  the  thennal 
speeds  of  hot  and  cold  electrons.  In  other  words  we  require  that  Vh  □  co/k  □  Vc .  This  means  that 
we  need  to  keep  C,c  □  1  and  <£h  D  1 .  The  plasma  dispersion  function  has  the  following  asymptotic 
expansion  for  |  C,  |  □  1: 

Z(C)D  -|{l+^  +  ^  +  ^  +  --j  (43) 

where 


■0  oft 

<T-jl  (44) 

2  r"  c _ l 

with  C,  =  C,'  +  i^"  .  In  the  other  limiting  case  when  |^|  □  1  we  use  the  following  power  series 
representation  for  Z  : 


=  - 2dl-3f  +  M  + 


(45) 


On  making  use  of  these  representations  in  (40)  and  equating  the  real  parts  we  have 


4  Resonant  damping  occurs  in  our  problem  due  to  the  energy  exchange  between  the  electron  acoustic  wave  and  particles  in  the  plasma  [T.H.  Stix,  {\it 
Waves  in  Plasma},  Springer- Verlag,  New  York,  1992].  Those  particles  having  velocities  less  than  the  phase  velocity  of  the  EAW  will  be  accelerated 
by  the  EAW.  As  commonly  done,  particle  velocities  are  modeled  as  a  Maxwellian  distribution  function.  Thus  there  are  more  particles  with  velocities 
less  than  the  wave  phase  velocity.  Consequently,  there  are  more  particles  gaining  energy  from  the  wave,  and  this  is  the  reason  for  plasma  waves  to  be 
damped. 
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1-*Re(*)  +  P7{1-'^Im(^)}  =  0  (46) 

Assuming  the  frequency  to  be  complex  by  letting  co  =  co'  +  ico" ,  and  imposing  the  condition  that 
\co"\U  (0  ,  we  obtain  the  following  solution 


CO  - 


CO,. 


l+k  n 


—2—2 
h  , 


(47) 


which  is  identical  to  the  result  that  we  obtained  using  the  hydrodynamic  approach.  We  follow  a 
similar  procedure  by  equating  the  imaginary  parts  of  (40)  and  obtain  the  following  expression  for 
damping  rate: 


CO  _ _ L 

co'  2  V  2  nh 


(48) 


We  thus  infer  that  the  ratio  njnh  has  to  be  very  small  in  order  to  keep  the  damping  rate 

acceptable.  However,  it  is  impractical  to  realize  this,  for  example,  in  reentry  applications.  Therefore, 
to  maintain  an  acceptable  level  of  EAW  energy  we  introduce  an  electron  beam  into  the  system. 

The  density  of  the  beam  electrons  can  be  small.  However  its  velocity  has  to  resonate  with  the  phase 
velocity  of  EAW.  With  the  addition  of  such  an  electron  beam  the  dispersion  equation  now  becomes 

1  +  itr{l  +  CMO)  +  P?{1  +  *Z(Q}  +  -£{1  +  &Z<6)}  =  0  (49) 

where 

^  =  (50) 


£,  =  4,-^  costf  (51) 

Vob  is  the  beam  velocity  and  6  is  the  orientation  angle  of  the  beam  with  respect  to  the  plasma 
sheath.  rh,  and  Vh  are  respectively  the  Debye  radius  and  thermal  speed  of  the  beam  electrons. 

We  keep  the  magnitude  of  £b  between  that  of  the  hot  and  cold  electrons.  It  is  not  possible  to  obtain 
a  closed-form  solution  to  the  above  dispersion  equation.  However,  one  can  get  an  estimate  of  the 
damping  characteristics  by  assuming  that  lO1  and  using  the  asymptotic  expansion.  Since 
nb  j nc  is  small  we  find  that  there  is  no  significant  change  in  the  real  part  of  the  frequency.  However, 
the  imaginary  part  of  the  frequency  is  different  and  is  given  as 


CO 

co' 


2  nu 


2  +  4^^ 

'5?  a  a 


-i 


(52) 


Notice  that  since  £b  is  much  smaller  than  £b  the  attenuation  in  the  scheme  with  an  electron  beam 
is  significantly  less  than  without  it. 


Indeed  these  estimates  are  based  on  asymptotic  relations.  In  general  one  needs  to  solve  the 
dispersion  equation  numerically  to  identify  the  regime  which  is  most  suitable  for  the  application  at 
hand.  Below  we  present  two  examples  to  illustrate  the  dispersion  characteristics  of  EAW.  In  the 

first  example  we  choose  the  following  set  of  parameters:  no  =  10Hcc_1 ,  nc  =  0.6no ,  nh  =  0.4«o , 
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Tc  =  1  eV ,  and  Th  =  20  eV  .  The  results  obtained  by  a  numerical  solution  of  (49)  are  shown  in 
Figure  4. 

The  complex  roots  of  the  dispersion  relation,  viz.,  the  real  and  imaginary  parts  of  the 
frequency,  correspond  to  the  real  frequency  obtained  by  our  hydrodynamic  approach  and  to  the 
wave  attenuation  due  to  resonant  damping.  In  Figure  4  we  have  shown  the  normalized  (normalized 
with  respect  to  cold  electron  plasma  frequency)  frequency  and  attenuation.  Note  that  the 
propagation  constant  is  normalized  using  the  Debye  radius  of  the  hot  electrons.  As  predicted  by  our 
analysis,  the  damping  is  very  high  even  though  we  have  chosen  the  parameters  favorable  to  the 
existence  of  EAW.  In  order  to  maintain  an  acceptable  level  of  EAW  we  have  introduced  into  the 
system  an  electron  beam  whose  velocity  resonates  with  the  phase  velocity  of  EAW.  To  illustrate  the 
significance  of  the  electron  beam  we  have  chosen  an  example  with  the  following  parameters: 

no  =10ncc-1,  nc  =0.9 na,  nh  =0.01  na,  nb  =0.09 n0,  Tc  =1  eV,  Th  =20  eV,  Th  =0.1  eV, 


Vh  =  6.5  x  108  cm/s .  The  solution  to  the  dispersion  equation  with  these  set  of  parameters  is  shown 

in  Figure  5.  Note  that  total  electron  density  here  is  the  same  as  in  Figure  4.  Yet  the  damping  in  this 
case  is  significantly  less  than  in  the  case  without  the  beam.  Since  our  interest  is  in  the  regime  when 
the  operating  frequency  is  much  less  than  the  cold  plasma  frequency,  i.e.  when  the  nonnalized 
frequency  ratio  is  very  small,  we  see  that  the  damping  is  fairly  small  here. 
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Figure  5.  Resonant  Damping  of  Electron  Acoustic  Waves  in  the  Presence  of  an  Electron  Beam. 
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3.5  Current  Sources  in  the  Two-Temperature  Electron  Population 


We  are  interested  in  studying  the  radiation  characteristics  of  a  source  distribution  embedded 
in  an  overdense  plasma  sheath.  If  it  is  sufficiently  thick  (much  larger  than  the  skin  depth)  there  will 
not  be  any  radiation  outside  the  plasma  sheath.  By  injecting  into  it  a  small  population  of  hot 
electrons  we  create  a  composite  medium  made  of  a  two-temperature  electron  population.  The 
embedded  antenna  can  excite  EAW  which  will  then  transform  into  EMW  at  the  boundary  of  the 
sheath  and  hence  radiate.  The  geometry  of  our  problem  is  shown  in  Figure  6.  The  lower  half-space 
is  composed  of  the  two-temperature  electron  population.  Embedded  in  this  medium  we  may  have  a 
distribution  of  current  sources.  We  are  interested  in  the  radiation  characteristics  of  these  current 
sources  in  free  space  (upper  half-space)  outside  the  plasma  medium.  To  facilitate  this  study  we 
derive  the  Green's  functions  for  the  problem.  We  denote  the  lower  half-space  as  Region  1  and  the 

upper  half-space  as  Region  0.  Thus  Gn  denotes  the  Green's  function  when  both  the  source  and 
observation  points  are  in  the  lower  half-space.  We  saw  that  in  our  two-temperature  electron 
population  there  can  exist  two  types  of  waves,  viz.,  EMW  and  EAW.  Hence  Gn  consists  of  two 
parts  as  given  below: 

Gjj  =  ciG  e  +  cxbW  V  Ga  (53) 

where  the  first  part  represents  EMW  and  the  second  part  represents  EAW.  More  specifically,  G 
represents  the  EMW  in  our  problem  and  is  governed  by 

V  x  V  x  Ge  -  k; Ge  =  IS(r  -  r')  (54) 

Ga  represents  the  EAW  for  our  problem  and  is  governed  by 


V  x  V  x  Ga  -  k~Ga  =  -S(r  -  r') 
The  constants  a  and  b  in  (53)  are  given  as 


where 


a  = 


e=  1—f- 

P  CO 


_  4  ni 


(55) 

(56) 

(57) 


a  is  a  constant  that  will  be  detennined  on  imposing  the  boundary  conditions.  The  propagation 
constants  ke  (for  EMW)  and  ka  (for  EAW)  are  given  as 


k  =k  |l-4 

co 


k„  = 


ylcO'-a 


(58) 


The  solution  for  our  problem  in  the  case  when  the  source  is  in  the  lower  half-space  and  the 
observation  point  in  the  upper  half-space  is  represented  by  the  Green's  function  G01 .  Noting  that  the 
upper  half-space  supports  only  EMW  the  governing  equation  is  given  as 

V  x  V  x  G01  -  k°G0l  -  0  (59) 
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where  k0  =  co/c  .  The  boundary  conditions  associated  with  these  Green's  functions  are  based  on  the 


following  requirements.  The  tangential  components  of  the  electric  and  magnetic  fields  must  be 
continuous  across  the  interface  z  =  0,  and  the  normal  components  of  the  velocity  vectors  of  the 
cold  and  hot  electrons  must  vanish  at  the  interface.  This  second  condition  ensures  that  the  charge 
distributions  of  the  two-temperature  electron  population  do  not  drift  away  and  are  confined  to  the 
geometrical  bounds  of  the  plasma  sheath.  These  requirements  lead  to  the  following  conditions  on 
our  Green's  functions 


(61) 

(62) 


(60) 


(63) 


where  r  and  rd  are  the  Debye  radii  of  the  cold  and  hot  electron  populations. 

A 

z 


Region  0 


free  space 


y 


Region  1 


2-T  electron  population 


Figure  6.  Geometry  of  the  Radiation  Problem.  The  antenna  is  located  on  the  z  -axis 
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3.6  Formulation  of  Green's  Functions 


First  we  construct  the  Green's  function  Gn  as  shown  in  (5 1).  As  given  there,  it  is  composed 
of  G  and  Ga  as  defined  in  (54)  and  (55).  The  solutions  for  G  and  Ga  can  be  readily  written  for 
our  geometry.  We  use  the  following  spectral  representations: 


G^(r,r')  =  ^Jc/k±i{(hV‘c-r  +  7?hVk«r)h+  +  (vVk*r  +  SvVk* r)  v+| 
G>,r')  =  ^  J  dk±  f(eik-T  + 


e~iKr'  (64) 


(65) 


where  h  and  v  are  unit  vectors  representing  horizontal  and  vertical  polarizations  of  EMW.  The 
superscripts  +  and  -  indicate  whether  the  corresponding  wave  is  propagating  upwards  or 
downwards.  R  and  S  are  the  reflection  coefficients  of  horizontally-  and  vertically-polarized  waves, 
respectively,  from  the  interface  z  =  0  when  the  EM  wave  is  incident  from  below.  ke  is  the  wave 
number  vector  of  the  electromagnetic  waves,  defined  as 

k*  =  k±  ±  ziri  U  =  ^k2p  +k“  (66) 

Notice  that  I]  is  positive  and  hence  the  z  component  of  ke  is  imaginary.  Thus  the  EMW  in  the 

lower  half-space  are  evanescent  in  our  frequency  regime  (overdense  case).  In  contrast,  the  wave 
number  vector  of  the  EAW  is  given  as 

k;=k  ±±zq,  ?„=A2-ki  (67) 

where  the  z  component  of  the  wave  number  vector  is  real  and  hence  corresponds  to  propagating 
waves.  The  superscript  >  on  G  and  Ga  indicates  that  the  representations  are  for  the  case  when 

z  >  z  .  The  corresponding  representations  for  the  case  when  z  <z  may  be  obtained  using  the 
following  symmetry  relations: 

Gj(r,r')  =  G^(r',r)  (68) 

G;(r,r')  =  G;(r',r)  (69) 

where  the  superscript  t  denotes  transpose.  Using  these  Green’s  functions  we  may  construct 
G  n  and  Gj)  .  However,  notice  that  we  only  need  Gj)  for  imposing  the  boundary  conditions  at  the 

z  =  0  interface.  We  construct  G01  noting  that  only  EMW  can  exist  in  the  upper  half-space.  Thus  the 


spectral  representation  for  G01  is  given  as 


01 

The  subscript  o  is  used  to  denote  the  polarization  vectors  in  the  upper  half-space.  The  propagation 
vector  in  the  upper  half-space  is  defined  as 

k„  =  kx  ±  zkoz  koz  =  yjk2-  k*  (71) 
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On  substituting  the  above  representations  of  Green's  functions  in  the  boundary  conditions  one  may 
obtain  solutions  for  the  various  coefficients.  Instead  of  trying  solve  this  general  boundary  value 
problem  it  is  much  easy  to  obtain  the  solution  that  corresponds  to  a  specific  source  current 
configuration. 
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4.  Results  and  Discussion 


4.1  Radiation  Pattern  of  a  Hertzian  Dipole 


With  the  availability  of  the  Green's  functions  we  may  calculate  the  radiation  characteristics 
of  any  source  current  distribution  in  our  composite  medium.  Since  we  are  interested  in  the  far- 11  eld 
radiation  pattern  in  the  upper  half-space  we  need  the  asymptotic  form  of  G01(r,r')  when  krU  1 . 


Evaluating  (70)  using  the  method  of  steepest  descent  we  obtain  the  following  asymptotic  result: 


G01(r,r')D 


k  k  +  7v+v+h~)e“'k"'r'  -4^7v+  (k+  +  5k  A  e 

o  oz  c  r  7  \  o  o  )  p  qa  o\  a  a  ) 


-ikTr' 


(72) 


Note  that  all  quantities  in  the  above  expression  are  evaluated  at  6 ,  which  is  the  elevation  angle 
associated  with  the  observation  point  r  .  Using  this  representation  of  the  Green's  function  we  can 
calculate  the  far-zone  radiating  fields  corresponding  to  arbitrary  source  current  distributions 
embedded  in  our  composite  medium.  Suppose  J  is  the  source  current  distribution  in  our  composite 

medium.  Its  far-zone  radiation  fields  in  the  upper  half-space  are  given  as 

o 

E(r)=jG01(r,r')-J(rVr'  (73) 


On  substituting  (72)  in  (53)  we  obtain  the  following  expressions  for  6  and  (j)  components  of  the 


electric  fields: 


EM0 


•J(r>-^V-"'r/r' 

+  5k  j  )  •  J  (r')e~ik±'r  e~iq“:  dr' 


(74) 

(75) 


For  sources  not  too  close  to  the  interface  we  hence  notice  that 


EA  r) 


is  exponentially  small. 


Indeed,  if  the  sources  are  very  close  to  the  interface  the  radiation  will  not  be  affected  much  by  the 
overdense  plasma  and  hence  there  is  no  need  to  create  a  two-temperature  electron  population.  It  is 
only  in  the  case  when  the  current  sources  are  well  away  from  the  interface  in  the  overdense  plasma 
that  we  need  to  inject  a  hot  electron  population  as  suggested  in  this  report  to  facilitate  low- 
frequency  signal  transmission  across  the  plasma.  Thus,  for  any  source  distribution  which  is  not  too 
close  to  the  interface  the  radiated  fields  are  entirely  0  -polarized.  In  this  sense  the  plasma  sheath 
may  hence  be  regarded  as  a  polarization  filter  that  filters  out  (j)  -polarized  signals. 


More  detailed  characteristics  of  the  radiation  pattern  may  be  inferred  only  when  we  know 
about  the  various  coefficients  that  appear  in  the  Green's  functions.  As  mentioned  earlier  these 
coefficients  are  obtained  by  imposing  the  boundary  conditions.  For  an  arbitrary  source  distribution 
it  is  fairly  complicated  to  evaluate  the  coefficients.  We  hence  focus  attention  on  two  specific  cases. 
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4.2  Vertical  Hertzian  Dipole 


In  this  case  we  choose  the  source  current  to  be 

J(r)  =  zJoS(r±)S(z  +  d ) 


(76) 


where  JQ  is  its  amplitude.  Substituting  this  in  (74)  and  (75)  we  obtain  the  following  solution  for  the 
far-zone  scattered  fields: 


E(r)  =  Y±  (1  ■ -  Sjjy 


(77) 


Further,  employing  the  boundary  conditions  we  obtain  the  following  solutions  for  the  various 
coefficients: 


fv  —  lkp  sin2  0 

-L  ma  l +k;R;, 

„  _  kp  cos  d-g-vpkp  cos  0o 

_L  kp  cos  0-%+vpkp  cos  0o 

Y  2 i Vpkp  sin 0\kp  cos 0 
±  _L  iD\_kp  sin  0+2  at]  vp  cos  0o 


Z  =  -ian  vp=f 

D±  =  kp  cos  0  -  %  +  v pkp  cos  0o 


(78) 

(79) 

(80) 

(81) 

(82) 


Notice  that  the  problem  has  azimuthal  symmetry  and  hence  the  above  expressions  are  independent 
of  azimuthal  angle.  We  have  used  the  subscript  _L  to  indicate  that  these  solutions  are  for  the 
vertically-polarized  Hertzian  dipole. 

To  illustrate  the  radiation  characteristics  of  the  vertical  dipole,  we  have  chosen  the  following 
numerical  example.  The  electron  density  is  chosen  as  1012[cc_1] ,  5  %  of  which  constitute  the  hot 
electron  population.  The  temperatures  of  the  cold  and  hot  electrons  are  0.2  eV  and  20  eV, 
respectively.  The  signal  frequency  of  the  antenna  is  taken  as  1  GHz.  Since  the  plasma  frequency  for 
this  example  is  around  9  GHz,  we  are  operating  in  the  overdense  region  and  hence  we  need  a 
scheme  like  the  two-temperature  electron  population  to  transmit  signals  across  the  plasma  sheath. 

As  explained  in  the  report  the  EAW  can  carry  low-frequency  microwave  signals  across  the 
overdense  region.  They  get  transformed  into  EMW  at  the  plasma  free-space  interface  and  hence 
radiate.  As  mentioned  earlier,  the  (f)  -polarized  components  of  the  electric  fields  are  exponentially 
small  and  hence  only  the  G  -polarized  components  radiate.  A  polar  plot  of  the  radiation  pattern  is 
shown  in  Figure  7.  The  radiation  pattern  of  the  vertical  dipole  embedded  in  the  plasma  is  shown  in 
red  and  that  in  free  space  is  shown  in  blue.  We  see  that  the  radiation  is  skewed  at  an  oblique  angle. 
This  is  typical  for  our  problem.  The  skew  angle  may  change  depending  on  the  parameters,  but  the 
radiation  pattern  is  always  found  to  be  skewed  at  some  angle. 

It  is  of  interest  to  compare  the  radiated  fields  for  this  overdense  case  with  those  of  an 
underdense  case.  When  the  electron  density  is  sufficiently  small  such  that  0)p  <  co ,  we  have  the 
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underdense  situation.  The  electron  density  in  our  example  is  108[cc_l] ,  which  corresponds  to  the 
plasma  frequency  if  around  90  MHz.  This  is  the  underdense  situation.  From  the  dispersion  relations 
(58)  we  see  that  ke  now  is  real  and  ka  is  imaginary.  Thus  the  EAW  are  now  exponentially  small 
and  it  is  now  the  EMW  that  carry  the  microwave  signals.  The  radiated  fields  are  given  as 

E(r)  =  Gko  YJ0e~q‘d  sin6>  (83) 

where 


Y  = 

2  kp  cos  0o 

kQ  cos  6+kp  cos  0o 

9  = 

{kl-kl  sin2^ 

(84) 

The  normalized  magnitude  of  the  radiated  fields  for  this  case  is  shown  (green)  in  Figure  7  for 
comparison. 
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Figure  7.  Comparison  of  radiation  patterns  of  a  vertical  dipole  in  plasma. 

We  next  turn  our  attention  to  the  other  important  case  when  the  dipole  is  oriented  parallel  to  the 
interface. 
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4.3  Horizontal  Hertzian  Dipole 


We  now  choose  the  source  to  be 

J(r)  =  zJoS(r±)S(z  +  d )  (85) 

where  p  is  a  unit  vector  on  the  horizontal  plane  (that  indicates  the  orientation  of  the  dipole). 
Substituting  this  in  (74)  and  (75)  we  obtain  the  following  solution  for  the  far-zone  radiated  fields. 

E(r)  =  -0kX  ^-g-Ya(l  +  Sa)j0  cos«(  (86) 

where  <j)  is  the  azimuth  angle  that  the  plane  of  observation  makes  with  the  p-z  plane. 

Further,  employing  the  boundary  conditions  we  obtain  the  following  solutions  for  the  various 
coefficients: 


kp  cos  0-ko  cos  0C 


X  = 


kp  cos  0+ko  cos  0o 

2k p  cos  6 
kn  cos  6+kn  cos  0n 


~  _  'X^cosd 

^  2^v/T  sin2  0 


cos  Qo  ±  (Tqav2p  cos2  9o  +  4 kp  sin2  9  cos  6 


(87) 

(88) 

(89) 


A  = 


iqa  cos  #(cos  O-Vp  cos  0o ) +a,7 j  sin2  0 
iqa  cos  #(cos  O+Vp  cos  0o  ^-aq  sin2  0 

2  iVpkp  cos2  6 
L\]kp  -2  aqkQ  cos  0  cos  0O 


Dt  -  iqa  cos  f?(cos  9  +  vp  cos  9o  -  dq  sin2  9 

*  =  1  +  k2aR2d 


(90) 

(91) 

(92) 

(93) 


The  subscript  □  is  used  to  indicate  that  these  solutions  correspond  to  the  case  of  the  horizontal 
dipole.  Notice  that  in  this  case  there  are  two  solutions  for  aQ .  For  the  parameters  of  our  problem  it 
turns  out  that  the  solution  with  the  positive  value  of  the  discriminant  is  nonphysical.  Hence  we  take 
the  value  for  au  with  the  negative  discriminant  as  the  solution.  As  before,  the  (j)  -components  of  the 

electric  fields  are  exponentially  small;  only  their  9  -components  are  radiated.  Figure  8  shows  the 
radiation  pattern  in  the  vertical  plane  parallel  to  the  dipole.  As  before,  the  radiation  pattern  for  the 
horizontally-oriented  dipole  embedded  in  the  overdense  plasma  is  shown  in  red.  The  corresponding 
results  for  the  dipole  in  free  space  are  shown  in  blue.  Actually  the  blue  curve  is  not  visible  because 
it  is  overshadowed  by  the  green  curve  that  represents  the  underdense  case.  For  illustration,  we  have 
chosen  the  same  parameters  as  for  the  vertical  dipole.  We  see  that  there  is  a  significant  skew 

towards  9  =  90“ .  More  importantly  there  is  a  complete  absence  of  radiation  in  the  broad-side 
region.  In  other  words,  the  radiation  is  close  to  the  end-fire  direction.  Thus  this  configuration  is  not 
suitable  for  constructing  phased  arrays  intended  for  beam  sweeping  around  the  broad  side. 
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Figure  8.  Comparison  of  radiation  pattern  of  horizontal  dipole  in  plasma. 

We  now  look  at  the  underdense  case  for  the  horizontal  dipole.  Again  the  EAW  are  exponentially 
small  and  hence  our  microwave  signals  are  carried  entirely  by  EMW.  The  radiated  field  in  this  case 
is  given  as 


E(r)  =  kakoz  J0  \~9Y cos0cos^  +  ^Xsin^}  (94) 

where  X  and  Y  are  given  earlier.  As  expected,  the  radiated  fields  contain  both  0  -  and  (f)  -polarized 


waves.  The  normalized  amplitude  of  the  electric  fields  in  the  p-z  plane  is  shown  (green)  in 
Figure  8  for  comparison.  Note  that  the  result  coincides  with  the  normalized  amplitude  of  the  electric 
field  for  the  dipole  in  free-space. 

In  the  next  section  we  turn  our  attention  to  the  radiation  characteristics  of  an  arbitrarily-oriented 
Hertzian  dipole. 
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4.4  Hertzian  Dipole  with  Arbitrary  Orientation 


The  source  current  of  the  Hertzian  dipole  with  arbitrary  orientation  is  represented  as 

3{r)  =  [pJp  +  zJz')S{rA)8{z  +  d)  (95) 

where  J p  and  J  are  the  components  of  the  source  current.  As  before,  the  amplitude  of  this  current 

is  J 0  and  its  direction  is  denoted  by  the  unit  vector  J  .  Using  the  results  for  the  horizontally-  and 
vertically-oriented  dipoles  we  may  construct  the  Green's  function  for  our  problem  as 

G01(r,r')  =  GL-zz  +  GQ-  pp  (96) 


where 


G,(r,r')  =  iirj^kJ 


-ik  r 


tl-w 


+ 


)e-v--‘+aYpf  v:(k:+S,k;)e-*:''}  (97) 


and  p  =  { _L,  □}  .  Thus  the  far-zone  radiated  field  from  an  arbitrarily-oriented  Hertzian  dipole 
located  at  2  =  —d  is  given  as 

E(r)  =  G01(r;(W)-J  (98) 

Thus  we  can  calculate  the  radiation  characteristics  of  an  arbitrary  source  distribution  embedded  in 
the  overdense  plasma. 


We  hence  see  that  a  two-temperature  electron  population  may  be  used  to  transmit  low-frequency 
signals  across  an  overdense  plasma.  However,  there  are  peculiarities  and  limitations.  We  found  that 
only  6  -polarized  waves  are  radiated.  Further  the  radiation  is  skewed  towards  large  6  values.  This 
is  not  a  major  set-back  for  the  vertically-oriented  dipole.  However,  the  horizontally-oriented  dipole 
has  serious  limitations  in  its  ability  to  be  used  in  phased  arrays  intended  for  scanning  in  the 
broadside. 
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5.  Conclusion 


We  have  studied  in  this  report  the  radiation  characteristics  of  a  source  distribution  in  an  overdense 
plasma  composed  of  a  two-temperature  population.  Indeed,  low-frequency  EM  signals  will  be 
completely  cut  off  in  this  case  and  hence  will  not  be  able  to  radiate.  However,  the  two-temperature 
electron  population  supports  EAW  whose  dispersion  relation  complements  that  of  the  EMW  in  the 
sense  that  EAW  propagate  in  the  frequency  domain  where  EMW  are  cut  off  and  vice  versa.  Thus  in 
principle  one  may  employ  EAW  to  carry  low-frequency  microwave  signals  across  an  overdense 
plasma.  Unfortunately  resonant  damping  can  seriously  limit  this  scheme.  However,  we  find  that 
with  the  introduction  of  an  electron  beam  into  the  system  it  is  possible  to  keep  the  damping  rate  to 
an  acceptable  level  for  our  application5. 

We  have  ignored  collisions  in  our  analysis  primarily  to  focus  on  the  concepts  of  EAW  without  the 
complexities  associated  with  the  inclusion  of  collisions.  We  admit  that  there  will  be  collisions 
within  the  plasma  sheath.  However,  we  find  that  in  the  domain  of  our  interest  the  collision 
frequency  is  much  smaller  than  the  plasma  frequency  [4].  As  with  any  plasma  wave  process  the 
electron  acoustic  waves  will  undergo  certain  amount  of  attenuation  because  of  collisions.  However 
attenuation  due  to  resonant  damping  is  far  more  significant  than  that  due  to  collisions.  Even  in  the 
region  of  resonant  damping  our  analysis  in  Section  V  shows  that  the  real  part  of  frequency  of  the 
EAW  is  unaffected  (see  eqn  (47)).  The  primary  consequence  of  resonant  damping  is  increased 
attenuation.  Therefore,  collisions  will  not  alter  the  fundamental  character  of  EAW.  In  particular,  the 
dispersion  relation  of  EAW  as  shown  in  Figure  2  will  remain  unaffected.  Our  initial  study  confirms 
this.  Collisions  will  only  introduce  additional  attenuation  of  the  electron  acoustic  waves. 


To  facilitate  our  study  of  the  radiation  characteristics  of  an  arbitrary  source  distribution  embedded 
in  this  composite  medium  we  derived  the  Green's  function  for  our  problem.  We  hence  find  that  only 
6  -polarized  waves  can  radiate  from  our  system.  To  understand  the  characteristics  of  our  system  we 
computed  the  radiation  patterns  of  vertically-  and  horizontally-oriented  dipoles  embedded  in  our 
medium.  It  is  found  that  the  radiation  pattern  of  horizontally-polarized  waves  has  two  lobes  (in 
contrast  to  the  underdense  case)  which  are  strongly  skewed  towards  the  horizon.  This  property 
makes  this  configuration  unsuitable  for  constructing  phased  arrays  intended  for  beam  sweeping 
around  the  broad  side.  To  summarize,  we  have  found  that  one  can  use  EAW  to  transmit  low- 
frequency  signals  across  an  overdense  plasma.  However,  there  are  peculiarities  and  limitations  that 
one  should  be  aware  of  before  embarking  on  the  design  of  sensor  systems  aboard  hypersonic 
vehicles.  For  instance,  we  have  assumed  that  we  can  create  a  stable  two-temperature  electron 
population  by  injecting  a  small  percentage  of  hot  electrons  into  the  cold  electron  population  that 
constitutes  the  overdense  hypersonic  flow  field.  The  hot  electrons  may  be  injected  by  means  of  an 
electron  gun  on  board  the  hypersonic  vehicle  like  it  was  done  in  the  CHARGE-2B  experiments 
[28].  However,  more  theoretical  and  experimental  studies  need  to  be  carried  out  on  how  to 


5  Indeed  there  are  engineering  issues  that  need  to  be  addressed  while  implementing  the  ideas  of  this  scheme  in  practice 
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successfully  implement  this  scheme  for  our  application.  Another  issue  concerns  the  transformation 
of  EAW  to  EMW.  We  are  aware  that  this  process  is  not  as  efficient  as  we  would  like  it  to  be.  There 
are  limitations  on  the  power  of  the  antennas  below  which  one  needs  to  operate  to  implement  this 
scheme.  Investigation  of  these  issues  and  finding  ways  to  mitigate  them  requires  further  detailed 
study  far  beyond  the  scope  of  this  report.  We  continue  to  actively  work  on  these  problems  details  of 
which  will  be  reported  in  future  publications. 
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List  of  Acronyms,  Abbreviations,  and  Symbols 


Acronym  Description 


EMW 

Electromagnetic  waves 

EAW 

Electron  Acoustic  Waves 

E 

Electric  field 

J 

Source  current 

no 

Electron  number  density 

K 

Wave  number  of  EMW 

K 

Wave  number  of  EAW 

% 

Plasma  Frequency 

z 

Plasma  Dispersion  Function 

rh’rc 

Debye  radii  of  hot  and  cold  electrons 
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